REPORT  DOCUMENTATION  PAGE 


fcm 

OM»  NO.  0704-0t99 


1.  A6£NCY  USE  ONLY  (Uiv*  titank)  1 2.  REMRT  OATE 

3-31-96 


4.  mu  AND  SUBTITU  , ,  . 

FINAL  REPORT  "Ndn-LineacrOptical  Effects  in 
Molecules"  AF0SR-F4962O-93-1-O1' 18 

6.  AUTHORS) 

Rodney  J.  Bartlett 

7.  KiVORMING  ORGANIZATION  NAME(S)  AND  AOOR£SS(ES) 

University  of  Florida 
219  Grinter  Hall 
Gainesville,  FL  32611 


3.  REPORT  TYPE  AND  OATES  COVERED 

.FINAL  REPORT  (2-1-93/1-31-96) 


is.  aiNOtNG  NUMBERS 


aso'i  /  !=■' 


AFOSR-TR-96 


12b.  DISTRIBUTION  CODE 


12a.  DISTRIBUTION /AVAILABiUTY  STATEMENT 

Approved  -for  public  release 
Distribution  is  unlimited. 


13.  ABSTRACT  (Msmmwn  2Q0viotdsf  analyse  the  evolution  of  fitst  principle  molecular  theory  m  obtaining 

reliable  values  of  molecular  hyperpolarizabilities.  Such  quantities  place 
severe  demands  on  the  capability  of  quantum  chemistry,  as  large  basis 
sets,  frequency  dependence,  and  high  levels  of  electron  correlation  are 
all  shown  to  be  essential  in  obtaining  observed,  gas  phase  hyperpolar- 
'  ‘  izabilities  to  within  an  error  of  10%.  Ab  initio  Hartree  Fock  results  are 
typically  in  error  by  nearly  a  factor  of  two,  while  the  errors  due  to  the 
neglect  of  frequency  dependence  average  ~  10-30%  depending  upon 
the  particular  process.  It  is  also  shown  that  for  small  molecules,  stan¬ 
dard  semi-empirical  approaches,  like  INDO  and  INDO/S,  will  often  not 
even  give  the  correct  sign  for  hyperpolarizabilities.  It  is  demonstrated 
that  using  state-of-the-art  correlated,  frequency  dependent  methods,  it 
"  "  is  possible  to  provide  results  to  within  10%.  Excluding  any  one  of  the 

essential  elements,  though,  destroys  the  agreement. 


15.  NUMBER  Of  RAGES 

9 


16.  PRICE  CODE 


20.  UMITATION  OF  ABSTRACT 

Unlimited 


Standard  Form  298  (Rev.  2-89) 

^rvicrtM  AMU  U4. 

m-loa 


14.  SUBJECT  TERMS 

-- 

- 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

Unclassified 

19.  SECURITY  CDVSSifiCATION 

OF  ABSTRACT 

Unclassified 

MSN  7540^1  •2aQ-S500 


3TXC  mSPlOTSD  I 


FINAL  REPORT 


Non-Linear  Optical  Effects  in  Molecules 
AFOSR-F49620-93-1-0118 


to  the  Air  Force  Office  of  Scientific  Research 
Chemistry  and  Materials  Science 
02/01/93  through  01/31/96 


Rodney  J.  Bartlett 
Quantum  Theory  Project 
University  of  Florida 
Gainesville,  Florida  32611-8435 


^  5  M/IH  isg§ 


TABLE  OF  CONTENTS 


Abstract 

2 

Can  Quantum  Chemistry  Provide  Reliable  Molecular 
Hyperpolarizabilities? 

3-36 

References 

NLO  Publications  (1979-present) 

36-37 

Invited  Presentations 

37-38 

University  CoUoquia  Presented 

39 

Postdoctoral  Researchers 

39 

Graduate  Students 

39 

Professional  Honors 

39-40 

1 


Can  Quantum  Chemistry  Provide 
Reliable  Molecular  Hyperpolarizabilities? 

Rodney  J.  Bartlett  and  Hideo  Sekino 
Quantum  Theory  Project 
Departments  of  Chemistry  and  Physics 
University  of  Florida 
Gainesville,  Florida,  USA  32611-8435 


In  the  late  70’ s,  one  of  us  (RIB)  became  interested  in  molecular  hyperpolarizabilities 
as  the  essential  element  in  non-linear  optics  (NLO)  when  Gordon  Wepfer  at  the  Air 
Force  Office  of  Scientific  Research  called  attention  to  J.  F.  Ward’s  results  from  dc- 
induced  second  harmonic  generation  (dcSHG)  experiments  and  the  enormity  of  their 
discrepancy  with  theoretical  results.  Tables  I  and  II  are  extracted  from  a  slightly  later, 
1979,  paper  of  Ward  and  Miller  (1)  demonstrating  the  problem. 

Not  only  was  the  existing  theory  results  of  the  time  typically  in  error  by  a  factor 
of  3  to  5  in  magnitude  for  the  electric  susceptibility,  but  even  the  signs  were 

frequently  wrong.  (The  sign  is  positive  if  the  measured  quantity  is  positive, 
where  //  is  the  permanent  dipole  moment.)  When  experimental  numbers  were  obtained 
by  other  techniques,  which  should  only  differ  by  different  dispersion  effects  that  are 
typically  less  than  10%;  they,  too,  had  little  correspondence  with  the  dcSHG  data. 
Note  from  Table  I  the  -3500  value  for  NH3  from  refractivity  virial  data  compared  to 
-209+5  for  the  dcSHG, 


.(2) 


Table  I  in  units  of  10'^^  esu/molecule.  Theoretical  values  from  various 
molecular  orbital  calculations  —  semi-empirical  (SE),  uncoupled  Hartree-Fock  (HF), 

and  coupled  Hartree-Fock  (CHF)  —  and  a  single  other  experimental  value  are  included 

( 2^ 

for  comparison.  The  sign  of  //Xy  is  unambiguously  determined  by  the  experiment 
and  is  independent  of  the  sense  chosen  for  the  molecular  z  axis.  This  table  extracted 
from  J.  F.  Ward  and  C.  K.  Miller  (7). 


THEORY 

Semi- 

Empirical 

Uncoupled 

HF 

Coupled 

HF 

x\\^  SHG 


Other 

EXP 


C-0+ 


0.112 

±0.005 


+129  ±14 


-43.5'’ 

(+95)*' 


+879*= 

+420*= 

+387*= 

-438® 


Table  I  in  units  of  10'^^  esu/molecule.  Theoretical  values  from  various 
molecular  orbital  calculations  —  semi-empirical  (SE),  uncoupled  Hartree-Fock  (HF), 

and  coupled  Hartree-Fock  (CHF)  —  and  a  single  other  experimental  value  are  included 

(2) 

for  comparison.  The  sign  of  //Xy  is  unambiguously  determined  by  the  experiment 
and  is  independent  of  the  sense  chosen  for  the  molecular  z  axis.  This  table  extracted 
from  J.  F.  Ward  and  C.  K.  Miller  (i). 


xf  SHG 

THEORY 

Other 

EXP 

Semi- 

Empirical 

Uncoupled 

HF 

Coupled 

HF 

C-0+ 

0.112 

+0.005 

+129  ±14 

+879‘= 

+420^= 

+387‘= 

-438® 

N-0+ 

0.15872 

±0.00002 

+147  ±17 

+47.7'’ 

H2+S- 

0.974 

+0.005 

-43  ±9 

N-H3+ 

1.474 

±0.009 

-209  ±5 

56.4’’ 

-44.4^ 

-65.1^ 

-19.0^ 

-15.6’’ 

-40.8’’ 

-35008 

H2^0- 

1.86 

±0.02 

-94  ±4 

120^ 

90.6 

-52.5^ 

-79.2^ 

-21.9’^ 

-51.6’’ 

-48.0’’ 

^  Electric  dipole  moments  in  Debye  units  from  Landolt-Bemstein,  Zahlenwerte 
und  Funktionen,  Neue  Serie,  Vols.  II/4  and  II/6  (Springer- Verlag,  Berlin)  and  Ref. 
18.  The  sense  of  the  NO  moment  is  suggested  by  F.  P.  Billingsley  H,  J.  Chem.  Phys. 
63  2267  (1975);  62,  864  (1975). 

^  N.  S.  Hush  and  M.  L.  Williams,  Theoret.  Chim.  Acta  (Berlin)  (25),  346 
(1972).  Signs  for  NH3  and  H2O  are  ambiguous. 

^  J.  M.  O’Hare  and  R.  P.  Hurst,  J.  Chem.  Phys.  {46),  2356  (1967). 

(0;0,0)  from  A.  D.  McLean  and  M.  Yoshimine,  J.  Chem.  Phys.  (46), 
3682  (1967). 

®  S.  P.  Liebmann  and  J.  W.  Moskowitz,  J.  Chem.  Phys.  54,  3622  (1971). 
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^  P.  Lazzeretti  and  R.  Zanasi,  Chem.  Phys.  Lett.  (59),  323  (1976)  Note  added 
in  proof.  Recent  reconsideration  of  the  sign  of  these  entries  yields  the  negative  signs 
now  shown  here  —  P.  Lazzeretti  (private  communication).  Overall  consistency  is 
substantially  improved  by  this  change. 

®  Refractivity  virial  data  from  A.  R.  Blythe,  J.  D.  Lambert,  P.  J.  Petter  and  H. 
Spoel,  Proc.  R.  Soc.  (London)  A  255,  427  (1960). 

^  G.  P.  Arrighini,  M.  Maestro  and  R.  Moccia,  Symp.  Farad.  Soc.  2,  48  (1968). 


In  Table  H,  the  discrepancy  among  different  experiments  is  even  more  apparent 
for  X||  ’  compared  to  other  dcSHG  values,  three  wave  mixing  (TWM),  Kerr  effect, 
and  third  harmonic  generation  (THG)  experiments.  The  theory  is  even  off  by  about  a 
factor  of  2  from  a  static  value  for  H2.  The  one  reasonable  theoretical  number  they  cite 
is  our  static,  correlated  value  for  N2,  one  of  our  first  (unpublished),  that  is  beginning 
to  be  in  reasonable  agreement  with 
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Table  II  in  units  of  lO”^^  esu/molecule.  Values  from  other  dc-electric 
field-induced  second-harmonic  generation  (dcSHG),  three-wave  mixing  (TWM), 
Kerr  effect,  and  third-harmonic  generation  (THG)  experiments,  along  with  theoretical 
results,  are  included  for  comparison.  This  table  extracted  from  J.  F.  Ward  and  C. 
K.  Miller  (7). 


y(3) 

X|| 

dcSHG 

Ward& 

Miller 

dcSH& 

TWM^ 

Kerr 

THG^= 

Theory 

H2 

65.2  ±0.8 

79 

- 

47^  ±5 

80  ±12 

34® 

N2 

86.6  ±1.0 

- 

104 

120^  ±10 

107  ±17 

718 

02 

95.3  ±1.6 

110 

100 

- 

- 

- 

CO2 

111.9  ±1.3 

- 

192 

750^  ±160 

156  ±23 

CO 

144  ±4 

- 

138 

- 

- 

- 

NO 

235  ±7 

- 

322 

- 

- 

- 

H2S 

865  ±22 

- 

- 

- 

- 

- 

NH3 

511  ±9 

- 

- 

- 

- 

- 

H2O 

194  ±10 

- 

- 

- 

- 

- 

®  Data  from  G.  Mayer,  C.  R.  Acad.  Sci.  B276,  54  (1968)  and  G.  Hauchecome,  F. 
Kerberve  and  G.  Mayer,  J.  Phys.  (Paris)  32,  47  (1971),  normalized  using  the  dcSHG 
coefficient  for  argon  from  R.  S.  Finn  and  J.  F.  Ward,  Phys.  Rev.  Lett.  26,  285  (1971). 

Data  from  W.  G.  Rado,  Phys.  Lett.  11,  123  (1967),  normalized  using  the 
dcSHG  coefficient  for  argon  from  R.  S.  Finn  and  J,  F.  Ward,  Phys.  Rev.  Lett.  26, 
285  (1971). 

^  J.  F.  Ward  and  G.  H.  C.  New,  Phys.  Rev.  185,  579  (1969). 

A.  D.  Buckingham  and  B.  J.  Orr,  Proc.  R.  Soc.  (London)  A305,  259  (1968). 

®  (0;0,0,0)  from  A.  D.  McLean  and  M.  Yoshimine,  J.  Chem.  Phys.  46, 

3682  (1967). 

^  A.  D.  Buckingham,  M.  P.  Bogard,  D.  A.  Dunmur,  C.  P.  Hobbs  and  B.  J.  Orr, 
Trans.  Fora,  Soc.  66,  1548  (1970). 

®  Xzzzz  (0;0,0,0)  from  R.  J.  Bartlett  and  G.  D.  Purvis  (private  communication). 

This  chapter,  which  is  intended  to  be  useful  to  experimentalists  who  are  trying 
to  assess  the  reliability  of  the  theory,  but  also  to  theoreticians,  as  it  provides  an 
overview  of  the  theory  which  can  and  has  been  used,  addresses  the  demands  that 
hyperpolarizabilities  place  on  first-principle  ab  initio  electronic  structure  theory;  and 
the  level  required  for  the  adequate  evaluation  of  molecular  hyperpolarizabilities.  After 
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doing  some  basic  elementary  perturbation  theory  that  underlies  all  that  is  done,  and 
which  provides  definitions  and  a  framework  for  discussion,  we  will  explain  the  various 
levels  of  quantum  chemical  application,  illustrated  by  numerical  results  we  have 
obtained  that  emphasize  the  various  approximations  employed  in  electronic  structure 
in  its  application  to  hyperpolarizabilities.  We  will  also  demonstrate  how  the  theory 
has  necessarily  evolved  to  provide  a  realistic  treatment  of  such  properties. 

The  order  of  presentation  will  follow  our  own  odyssey  that  led  us  to  introduce 
electron  correlation  into  molecular  hyperpolarizabilities  (2,3),  using,  then  new,  many- 
body  perturbation  theory  (MBPT)  methods  «5);  to  make  a  prediction  for  the 
FH  molecule,  that  had  later  ramifications  (6,7);  to  use  new  coupled-cluster  (CC) 
correlated  methods  including  those  augmented  by  triples  (8,9);  to  explore  vibrational 
polarizabilities  in  static  electric  fields  (10);  and  to  introduce  frequency  dependence  by 
developing  analytical  derivative  time-dependent  Hartree  Fock  (TDHF)  theory  (11,12). 
The  latter,  which  enables  a  first  (decoupled)  treatment  of  both  the  essential  frequency 
dependent  and  correlation  aspects  of  the  problem  (13),  culminates  in  a  uniform 
study  of  ten  molecules  (14).  The  correlation  calculations  were  assisted  by  parallel 
developments  in  analytical  derivative  CC/MBPT  theory  (15-18).  Finally,  we  coupled 
correlation  and  frequency  dependence  in  the  new  equation-of-motion  (EOM)  CC 
method  for  hyperpolarizabilities  (19,20).  By  taking  this  evolutionary  viewpoint,  we 
hope  the  current  review  will  complement  the  other  recent  excellent  reviews  (21-23) 
on  the  topic. 

The  theory  also  offers  different,  conceptual  viewpoints  on  hyperpolarizability 
evaluation  and  interpretation,  and  this  degree  of  flexibility  should  be  better  appreciated 
when  comparing  theoretical  numbers  and  addressing  NLO  design  criteria.  We  will 
conclude  with  some  recommendations  for  some  future  developments  in  the  continuing 
evaluation  of  predictive  quantum  chemical  methods  for  NLO  material  design. 

Perturbation  Theory  of  Molecular  Hyperpolarizabilities 

The  quantities  of  interest  are  the  electric  susceptibilities  and  In  the  gas  phase 
experiments  of  interest  here,  there  is  no  particular  distinction  between  macroscopic 
and  microscopic  susceptibilities  and  the  hyperpolarizabilities,  as  they  are  simply 
related.  Those  of  particular  interest  are  obtained  from  electric  field  (dc)  induced 
second  harmonic  generation  experiments  (1).  Namely,  a  sample  of  gas  is  subjected 
to  a  dc  field  (eo)  and  an  optical  electric  field  £u)(e'-^^  at  frequency  w,  to 

induce  a  dipole  moment  at  frequency  2uj.  Allowing  for  the  various  manipulations 
required  to  relate  the  molecular  quantities  to  the  laboratory  fixed  observables, 

-  (1) 

x\(-2uj;o,u,w)  =  x||^^(-2a;;o,a;,a;)  -|- 

X||  (-2a;;a;,w)  for  second  harmonic  generation  (SHG)  is  separated  from 

/o\ 

X||  ’(-2uj;o,u},uj),  which  is  dcSHG,  by  studying  as  a  function  of  T.  Numbers 
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relative  to  a  standard,  typically  He  gas,  are  obtained.  Revisions  of  the  He  reference 
value  cause  some  slight  rescaling  of  the  experimental  numbers  (22).  Furthermore, 
we  relate  the  susceptibilities  to  the  molecular  hyperpolarizabilities  via, 

= /3||(-2a;;w,a;)/2 

( s') 

X||  H-2a;;o,a;,a;)  =  7||(-2a;;o,w,a;)/6 

where  the  Einstein  summation  convention  is  employed,  meaning  all  repeated  indices 
are  summed,  i.e.  =  'dxx+''^yy  +  '^zz-  The  parallel  designation  (||)  means  measured 
parallel  to  the  dc  field,  while  the  (±)  component  can  be  similarly  obtained.  As 
and  7||  are  properties  of  molecules,  they  constitute  the  objective  for  first-principle 
quantum  mechanical  evaluation;  the  subject  of  this  chapter. 

The  basic  idea  underlying  any  treatment  of  a  molecular  response  to  an  electric 
field,  static  or  oscillatory,  is  the  solution  of  the  molecular  Schrodinger  equation  in  the 
presence  of  the  field.  We  will  first  consider  the  static  case,  generalizing  the  approach 

for  frequency  dependence  in  Section  7.  Our  perturbation  Eq-W  =  ein  is  given 

i 

by  the  interaction  of  an  electron  Si  at  position  fi  with  a  static  electric  field, 

£o  =  +  £yoj  +  £zok,  (3) 

where  £xo,  £yo,  and  Ezo  are  the  field  strengths  in  the  x,  y,  and  z  directions  of  magnitude 
|f:|.  In  atomic  units,  e  becomes  minus  unity,  so  we  have  a  Hamiltonian, 

n{£o)  ^Ho+So-W  =  Ho- £o  •  (4) 

i 

where  Ho  consists  of  the  usual  kinetic  energy,  —  potential  energy, 

i 

+  +  "fe/’  contposed  of  the  electron  nuclear  attraction,  the  two- 

electron,  electron-electron  repulsion  operator,  and  the  proton-proton  repulsion.  Greek 
letters  indicate  atoms  and  ij  electrons.  W  contains  the  negative  sign.  The  solution 
to  the  time  independent  Schrodinger  equation 

'Hork  =  EM  (5) 

provides  the  ground,  V’o  =  C  and  excited  states,  {ipl},  and  associated  eigenvalues, 
{E^},  for  the  unperturbed  molecule.  Notice  these  are  the  exact,  many-particle 
unperturbed  states. 

We  now  seek  a  solution  for  the  perturbed  Hamiltonian  in  its  ground  state 

=  E{£oM£o)  (6) 
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The  natural  way  to  attempt  a  solution  is  perturbation  theory.  Hence,  we  expand, 
H{£o),  and  E{eo)  in  a  perturbation  series  in  Eq.  For  the  time  being,  we  will 

assume  the  field  lies  solely  in  the  z  direction,  so  Exo  =  £yo  =  0  and  Ezo  =  £o-  Then, 


{Ho  +  EoW)  i^o  +  +  ■ .  •  j 

=  [e  +  EoE^^'^  +  £2£;(2)  +  . . .  j  +  £^^(1)  +  . . .  j 

As  the  equality  has  to  be  true  for  any  power  of  Eo  (i.e.  all  quantities  are  linearly 
independent),  gathering  terms  of  a  given  power  of  Eo  together,  we  obtain  for  e\\ 

{Eo  -  ={w-  (8) 

for  eI 

{Eo  -  no)i’^^^  =  (w-  -  E^^^o  (9) 

for  eI 

{Eo  -  no)i>^^^  =  (w-  -  E^^^o  -  (10) 

and  for  e^ 

{Eo  -  no)i)^^'^  =  (w-  ^(1)) -  £;(2)^(2)  (i  i) 

These  are  sufficient  to  take  us  through  the  electronic  dipole  moment,  //,  the  dipole 
polarizability  tensor,  a,  and  the  ^  and  7,  hyperpolarizability  tensors.  Multiplying  on 
the  left  by  (V'o|»  and  using  the  fact  that  {'tf)o\{Eo  —  T^o)  =  0,  we  have 

£;(i)  =  {^o\W\iPo)  (12) 


Similarly,  we  obtain 

£;(3)  =  {'tPo\w\i)^‘^^)  =  {i>^^^\w  -  (13) 

£;(4)  =  {'tPo\W\'lP^^'^)  =  {'lp^‘^^\Eo  -  -  £;(2)^^(i)|^(i)^ 

where  we  have  used  the  intermediate  normalization  condition,  =  0  for  m  ^ 

0.  The  second  form  comes  from  manipulations  that  demonstrate  the  2n+l  rule  that  an 
n*  order  wavefunction  will  determine  and  the  analogous  2n  rule  for  even  orders. 

Instead  of  straightforward  perturbation  theory,  we  can  also  derive  these  formulae 
from  explicit  differentiation  of  the  expectation  value  of  ’^(^o)  with  respect  to  Eo. 
From  equation  6,  we  have  the  expectation  value, 

£;(£o)(^(eo)|^(eo))  =  ('®'(£o)|W(£o)|^(£o))  (14) 
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and  differentiation  gives. 


{n{£o)  -  E{£o)) 


dso/ 


'd^ 

Mo 


dH  dE 
dso  deo 


(H{eo)  -  E{eo)) 


(15) 


As  we  want  to  evaluate  this  at  Eq  =  0,  using  {'Ho  —  Eo)\'ipo)  =  0,  and  its  complex 
conjugate  (c.c.),  gives 


dSo 


V’o 


J  eo=0 


m 

deo 


'>l>0 


=  {A\W\  i>o) 

£0=0  ' 


(16) 


This  is  simply  a  statement  of  the  simple  Hellman-Feynmann  theorem  that  the  deriva¬ 
tive  of  the  energy  is  given  as  the  expectation  value  of 
If  we  proceed  to  the  higher  derivatives,  we  will  obtain 

^'^E^  /  ,  I  ^  ^ 

dso  dso  I  dso 


^^Oj£o=0 


=  (V’o 


+C.C.  = 


m 

dso 


d'4>\ 


(17) 


or 


id'^E 
2  de^ 


/.  1 

dH 

II 

o 

II 

o 

to 

- 1 

deo 

d'lj) ' 
dso, 


(18) 


In  third  and  fourth  order. 


1  d^E 

^^J£„=0 


=  (i’o 


dH 

deo 


1  d'^E 

4! 


1  d^ 
2\del 

=  l'i>o 


d'lp 

deo 


dH  dE 
deo  deo 


d'lp 

deo 


dH 

deo 


1  d^i} 
3!  de^ 


(19) 


Clearly,  this  is  simply  perturbation  theory,  since  up  to  a  numerical  factor, 
~  ~  ~  5^’  words,  we  use  a 

McLaurin’s  series  expansion  for  £,  E{eo)  =  Eo+§^£o+^^^£l  +  ^^sl  +  . . .  and 
’i',  instead  of  a  straight  perturbation  series,  which  introduces  the  numerical  factors 
that  relate  the  perturbative  energies  to  the  derivatives  E^'^^  = 

We  now  define  the  electric  properties  as  the  electric  dipole  in  the  z  direction: 


dE 


=  -i^z  =  (i’o 


deo 

the  dipole  polarizability; 
d^E 


dH 

deo 


i’o 


i’o 


w 


i’o)=E^^^ 


del 


=  -OLzz  =  (  i’o 


deo 


di> 

deo 


^  =2 


the  first  hyperpolarizability; 


d^E 

del 


Pzzz  —  ^('4^0 


dH 

deo 


d^ 

del 


Mi’o 


(20) 


(21) 


)  =  3\E^^^  (22) 
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and  the  second  hyperpolarizability; 


^  =  -lzzzz  =  ^(^^o  W  ^(3)^  ^4,^(4)  (23) 


That  is, 


E(£o)  -  Eo  flzSzo  -  -^Olzz^lo  -  -^l^zzz^lo  “  -^Jzzzz^to  +  •  •  • 


which  gives  the  well  known  series. 

More  generally ,  when  all  components  are  considered,  with  the  Einstein  summation 
convention,  we  have 

E  =  E  -  HiEio  -  -^(Xijeioejo  -  -^Pijk^ioSjoSko  -  -^Jijki^io^jo^ko^io  -  •  •  •  (25) 

or  for  the  induced  dipole  moment, 

w(eo)  =  /ii(o)  +  aij£jo  +  ^PijkSjoSko  +  ^JijklSjoSkoSlo  +  ■■■  (26) 

Whether  the  numerical  factors  are  included  determines  the  choice  of  conventions  for 
polarizabilities.  Ward  and  co-workers  use  the  perturbative  definition,  meaning  that 
E  =  a,  £;(3)  =  and  ^(4)  =  7  ,  e.g.,  while  the  power  series  choice  directly 
associates  the  derivatives  with  the  polarizabilities,  obviating  the  numerical  factors. 
Since  this  choice  has  direct  correspondence  with  the  energy  derivatives,  it  appeals 
more  to  theoreticians. 

To  evaluate  the  polarizabilities,  we  require  a  knowledge  of  the  perturbed  wave- 
functions  (wavefunction  derivatives).  It  is  convenient  to  introduce  the  resolvent  op¬ 
erator, 

Ro  =  {Eo  -  Ho)~^Q  (27) 

where  Q  represents  the  projector  (Q^=Q)  of  all  functions  orthogonal  to  the  unper¬ 
turbed  reference,  ‘ipQ.  One  such  set  consists  of  all  other  eigenfunctions  of  Ho  (dis¬ 
crete  and  continuous),  making  <5  =  E  I  '^k^)Uk^  •  For  this  particular  set, 

k^o  *  '  ' 


{Eo  E^.  \'fpl)Q  —  Q[Eo  -  {^k^)-  By  virtue  of  excluding 

V’o.  the  inverse  Ro  operator  is  well  defined.  In  terms  of  the  resolvent, 

=  R^(w-  E^'^^^'iPo  =  RoWi>o 

_£;(!)  j^(l)  (28) 

etc.,  showing  that  all  order  wavefunctions  are  recursively  computed  from  the  prior 
ones.  Knowing  these  solutions,  we  can  compute  the  energy  from 

£;("+!)  =  {i;o\W\iP^^^)  (29) 
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or  explicitly,  in  a  few  cases. 


^(2)  =  {i,,\WRoW\i>o) 

£;(3)  =  {'^,\WRo{w  -  E^^'^^RoW\'il)o) 

=  {iPo\WRo{w  -  E^^^"jRo{w  -  -  E^^'^ii^olWRlWl'iPo) 

(30) 

etc. 

The  particular  choice  of  expansion  of  in  the  set  of  eigenfunctions  of  Ho  leads 
to  the  well  known  sum-over-state  (SOS)  formulas  for  the  z  components. 


-dz 


=E 

k^o 


l^z  =  {'lpo\z\i>o) 


Eo-El 


{o) 


-/3z 


Jzzzz  —  ^  ] 


E 

k.li^o 


(31a) 

(31b) 

(31c) 


(Eo  -  E^'^)  (Eo  -  e[^^) 

(Eo  -  4"^)  (Eo  -  4"^)  (K  - 
(V’o|^i4‘’^)(4‘’^|2|V'o) 


+o^ 


zz 


{eo  -  4'”)' 


(31d) 

The  other  tensor  elements  are  obtained  from  exchanging  x  and  y  with  z  in  all  distinct 
ways.  Henceforth,  we  will  recognize  that  all  practical  molecular  quantum  chemical 
calculations  employ  a  finite,  discrete  basis  set,  so  from  the  beginning  we  are  limited 
to  this  choice,  and  we  need  not  consider  any  continuum. 


Evaluation  of  Static  Hyperpolarizabilities 

From  the  above,  it  should  be  abundantly  clear  that  in  a  basis  the  derivative  approach, 
equations  20  to  23,  and  the  SOS  expressions,  equation  31,  are  simply  two  equivalent 
ways  of  expressing  polarizabilities.  Furthermore,  considering  that  the  choice  of 
excited  eigenfunctions  of  Ho  is  just  one  choice  for  a  complete  set  representation 
of  the  perturbed  wavefunctions,  tbe  above  SOS  forms  probably  attract  more 

significance  in  NLO  than  warranted.  For  example,  since  Eo  —  E^^  is  the  excitation 

energy,  and  since  (V’ol^lV'fc  /  is  the  z-component  of  the  transition  moment,  it  should 
be  possible  to  evaluate  a.^  from  equation  31b  purely  from  experiment  by  knowing 
all  electronic  excitation  energies  and  transition  moments  (including  those  for  the 
continuum  in  the  exact  case).  The  problem  is  knowing  them  all — a  very  large  number! 
Instead,  attempts  to  estimate  the  SOS  by  the  few  known  excitation  energies  and 


11 


transition  moments  is  likely  to  be  very  far  from  the  true  value  (see  our  contribution 
later  in  this  volume  (24)).  Note  for  the  ground  state,  all  contributions  to  from 
have  the  same  sign,  so  there  is  no  potential  cancellation  among  the  neglected 
contributions.  The  problem  is  further  compounded  for  P^zz  and  where  in  addition 
to  transition  moments  from  the  ground  to  excited  states,  it  is  necessary  to  know  the 
transition  moments  relating  two  excited  states,  and  that  information  is  hard  to  obtain. 
Note  that  can  have  either  sign.  77777  also  can  have  either  sign,  since  although 
the  lead  term  corresponds  to  {ipolEo  -  which  must  be  negative  (giving 

a  positive  contribution  to  'fzzzzX  the  second  term  is  positive,  attenuating  the  value 
of  77777.  Lacking  a  proof  that  the  magnitude  of  one  must  be  greater  than  the  other, 
either  sign  is  possible. 

From  the  above  we  have  two  viewpoints  on  the  evaluation  of  static  polarizabilities. 
We  can  either  evaluate  energy  derivatives  of  the  Schrodinger  equation  in  the  presence 
of  the  perturbation,  or  attempt  some  approximation  to  the  SOS.  Obviously,  the  former 
does  not  require  any  truncation.  (As  we  will  see  in  Section  8  later,  with  proper 
handling  neither  does  the  finite  basis  SOS.)  The  simplest  recipe  for  evaluation  of 
the  derivatives  is  to  use  what  is  called  the  finite-field  technique.  That  simply  means 
solve  the  Schrddinger  equation  in  the  presence  of  the  perturbation  by  choosing  a  small 
finite  value  for  So  of  0.001  a.u.,  e.g.  Adding  an  electric  field  quantity  to  Ho  gives  an 
unbounded  H(eo)  operator,  and  if  we  obtain  its  exact  solution,  the  lowest  energy  state 
would  be  the  field  ionized  state,  a  molecular  cation  plus  an  electron;  but  in  practice 
we  must  use  a  finite  basis  set  for  its  solution  which  is  effectively  like  putting  the 
molecule  into  a  box,  and  this  gives  a  valid  E{£o).  Repeating  the  procedure  at  £<,  = 
—0.001,  we  could  obtain  the  dipole  from  the  numerical  derivative. 


-dz  = 


d£o 


E{eo)  -  E(-£o) 

eo-»o  2£o 


(32) 


The  accuracy  depends  upon  the  size  of  the  finite  field  strength.  If  it  is  too  large,  the 
numerical  derivative  is  not  very  accurate,  while  if  too  small,  there  is  not  a  numerically 
significant  change  in  E(£o). 

To  obtain  all  p,  a,  ^  and  7,  we  need  several  more  points,  so  we  obtain 
expressions  like 


—  2  2£jo)j  +  —  £'(— £io)j  + 

7un4  =  ^\E{£io)  +  E{-£io)^  -  ^E{2£io)  -  E{-2£io)^  -  6E{o)  -f-  C>(4) 

See  (2)  for  others.  Note  each  expression  is  accurate  to  the  next  odd  (even)  order 
since  only  a,  7,  e  and  p,  S  are  interrelated  (2).  This  exclusion  of  the  next 
higher-order  contribution  greatly  helps  the  precision  of  such  a  calculation.  However, 
note  that  the  energy  needs  to  be  accurate  to  a  couple  of  significant  digits  better  than 
if  we  are  to  get  /0, That  is,  if  £jo  =0.001,  we  require  energies  to  be  10'^^  If 
£io  =0.01,  we  would  require  at  least  10'*.  For  7„„,  we  would  need  10'^°  to  lO'^"^.  As 
molecular  integrals  in  quantum  chemical  calculations  are  seldom  much  more  accurate 


12 


than  10'*^,  not  to  mention  other  parts  of  the  calculation,  finite  field  procedures  for 
hyperpolarizabilities  can  raise  serious  precision  problems. 

Another  problem  lies  in  the  proliferation  of  tensor  elements  in  ^  and  7.  Many 
energy  calculations  involving  field  strengths  in  different  directions  are  required  to 
evaluate  all  the  numerical  derivatives,  and  at  higher  levels  of  sophistication  these  are 
quite  expensive  calculations. 

The  solution  to  the  above  problem  is  to  analytically  evaluate  the  derivatives.  The 
simplest  is  —fiz  =  (V’ol^lV’o),  which  is  just  an  expectation  value.  For  the  others, 
analytical  means  that  while  solving  the  Schrddinger  equation  for  E,  we  also  directly 
obtain  all  components  of  the  derivatives,  ^  i  etc.  in  about  an 

equivalent  amount  of  time.  This  means  we  differentiate  before  evaluation  by  using 
equations  20  to  23  and  the  explicit  solutions  for  the  wavefunctions  and  derivatives, 
which  are  proportional  to  the  perturbed  wavefunctions  given  in  equation  28.  In  a 
substantial  formal  and  computational  achievement  of  30  years  duration,  primarily 
fueled  by  the  necessity  of  analytical  gradients  for  atomic  displacements  in  molecules 
{15,25,26),  such  analytical  procedures  have  been  developed  in  quantum  chemistry. 
Their  limitation  is  that  they  have  not  been  implemented  for  all  methods.  For  example, 
any-order  analytical  higher  derivatives  with  respect  to  electric  field  perturbations  have 
been  developed  for  the  Hartree-Fock  treatment  of  hyperpolarizabilities  by  exploiting 
the  recursive  nature  of  perturbation  theory  {11,27),  equations  28,  29.  For  correlated 
methods,  analytical  second-order  perturbative  theory,  [MBPT(2)],  derivatives  are 
available  for  a  and  ^  (75).  For  other  methods,  even  including  highly  sophisticated 
correlated  methods  like  coupled-cluster  theory  {28),  the  induced  dipole  /x(£o)  can  be 
evaluated  analytically,  from  which  numerical  derivatives  provide  a,  and  7  (13). 
In  this  way,  at  least  one  Eq  or  two  to  three  orders  of  magnitude  is  gained  in  the 
precision  of  and  7. 

Note  that  equation  31  represents  analytical  expressions,  too,  since  no  finite  field 
is  involved  in  their  evaluation.  The  latter  viewpoint  leads  to  the  analytical  evaluation 
of  the  (dynamic)  polarizability  using  the  equation-of-motion  (EOM)  CC  method 
{20).  Obviously,  it  would  be  ideal  to  be  able  to  analytically  evaluate  third  and 
fourth  derivatives  using  such  powerful  CC  correlated  methods,  but  the  theory  and 
implementation  has  not  yet  been  developed. 


Basis  Sets  and  Hartree-Fock  Theory 

Now  that  we  know  a  way  to  calculate  a  static  hyperpolarizability,  we  can  consider 
other  aspects  of  the  calculation.  The  first  approximation  to  consider  is  Hartree-Fock 
(HF),  self-consistent  field  (SCF)  theory.  That  is,  we  evaluate  the  energy  and  its 
derivatives  for  the  perturbed  Hamiltonian  by  obtaining  the  energetically  best  (lowest) 

single  determinant  solution,  =  A{(pi{l),  ^2(2)  •  •  •  (pnin)),  to  approximate  •0o.  In 

N 

the  absence  of  EoW,  that  means  FJrf  =  and  further,  Ho^Ho=  X)  /(O 

k=i 
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where  f  =  h  +  and  =  ^  /  y>j(2)^  ^t^y)j(2)dr.  The  effective  one- 

i=i 

particle  operator  is  an  average  over  the  two-particle  part.  The  orbitals  {pj} 

are  the  solutions  f{l)(pj{l)  =  ej(pj(l)  where  ej  is  the  HF-SCF  orbital  energy  (not 
electric  field  e  or  e  hyperpolarizability!).  Self-consistency  is  required  by  ensuring 
that  the  orbitals  {ipj}  used  in  in/are  self-consistent  with  the  solutions.  This 

model  provides  the  usual  molecular  orbital  (MO)  approximation  that  underlies  much 
of  our  conceptual  understanding  of  molecules.  This  averaging  procedure  introduces 
the  correlation  error,  which  pertains  to  electrons’  instantaneous  interaction  that  keeps 
them  apart.  It  corresponds  to  the  perturbation  which 

ifj  i 

we  will  consider  later  in  Section  6.  The  full  Hamiltonian  of  equation  41  is  thus 
n  =  Ho +  eoW  +  V. 

The  simplest  way  to  consider  doing  a  HF  calculation  of  a  polarizability  is 
the  finite-field  procedure.  That  means  that  we  compute  the  HF  approximation  to 
H  =  Ho  +  V  +  £oW  where  W  is  a  one-electron  operator  for  a  small  value  of  £o. 
Then,  we  will  obtain 


Eef{£o)  =  ($o(eo)|W(£:o)|$o(£o))  (34) 

where  $o(£^o)  =  -M}P£\{^)^e2{‘^)  •  •  •  ¥’en(^))  is  the  HF  wavefunction  and  its  compo¬ 
nent  orbitals  are  all  dependent  on  £o.  Furthermore, 

=  ej{£o)^jeil),  (35a) 

fe{l)  =  foil)  +  SoWil)  +  vlffil)  (35b) 

vfHl)  =  E/ <P%i2)^^:^cpj,i2)dT2  (35c) 


Obtaining  Eef{£o)  at  various  values  of  £o  will  provide  the  perturbed  energies  from 
which  the  numerical  derivatives  may  be  obtained.  This  finite-field  procedure  is 
frequently  called  “coupled  Hartree-Fock”  (CHF)  (29). 

Although  it  does  not  change  the  conceptual  content  or  the  numerical  values  (if 
done  carefully!),  the  much  more  computationally  convenient  analytical  equivalent, 
called  coupled  perturbed  Hartree-Fock  (CPHF),  can  be  developed  by  taking  the 
derivatives  before  evaluation  by  expanding  all  the  equations  in  perturbation  theory 

and  explicitly  solving  them  for  perturbed  orbitals,  pje  =  H - ,  orbital 

energies,  ej  =  ejo+£o^P  H - ,  and  using  -f  £oV^^^^^^  £lv^H if)  ^ - 

from  which  Eiif{£o)  =  +  ^ojE'hf  +  •  •  •  can  be  obtained. 

To  avoid  too  much  of  a  digression,  we  will  not  present  those  equations  here. 
Excellent  treatments  for  the  time-independent  case  are  given  elsewhere  (27,30,31). 
This  is  also  a  special  case  for  the  time-dependent,  TDHF,  approach  (11,32)  discussed 
in  Section  7.  Suffice  it  to  say  that  CPHF  calculations  are  preferable  to  CHF,  and 
several  implementations  are  available. 
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Even  at  the  HF  level,  though,  we  have  to  pay  close  attention  to  the  choice  of 
basis  set,  that  is  the  (usually)  contracted  Gaussian  atomic  orbital  (AO)  basis  used 
to  express  the  MO’s.  It  is  apparent  that  it  is  essential  to  have  a  large,  flexible 
and  polarized  basis  set  for  hyperpolarizabilities.  The  basis  must  correctly  describe 
matrix  elements  of  the  long-range  operator  r  (i.e.  x,y,z)  while  the  usual  AO  basis 
functions  have  been  selected  predominantly  to  describe  the  energy,  which  depends 
on  shorter-range  operators  like  ^  and  the  kinetic  energy  operator.  A  basis  like  double 
zeta  (DZ)  would  use  4s  and  two  sets  of  p  functions  to  describe  the  energy  of  a 
B,C,N,0,F,  atom  in  a  molecule,  2s  for  H.  To  provide  adequate  polarization  to  this 
basis,  particularly  to  describe  the  more  diffuse  and  directional  part  of  the  charge 
density,  we  need  at  least  one  or  two  sets  of  d-functions  and  probably  more  s  and  p 
functions.  The  POL+  basis  (14,33),  optimized  to  describe  polarizabilities  has  a  5s3p2d 
distribution  for  B-F,  and  3s2pld  for  H,  and  would  be  the  minimum  recommended  for 
most  hyperpolarizability  determinations.  The  influence  of  basis  set  is  illustrated  in 
the  behavior  of  the  components  of  a  and  7  for  ethylene,  shown  in  Table  HI. 

Table  HI  Comparison  of  Ab  Initio  SCF  Hyperpolarizabilities  as  a  Function  of  Basis 
Set  with  Semi-empirical  INDO  and  INDO/S  for  Ethylene  (a.u.). 


COMPONENT 

SEMI- 

EMPIRICAL 

AB  INITIO  SCF 

EXP 

INDO 

INDO/S 

STO-3G 

DZ 

63I+PD 

POL+ 

Otxx 

19.9 

31.7 

11.45 

33.6 

36.0 

36.4 

ayy 

15.9 

18.1 

0.75 

18.0 

22.9 

24.6 

2.8 

3.7 

2.84 

8.6 

19.4 

23.1 

a 

12.9 

17.8 

5.01 

18.4 

26.1 

28.0 

28.7 

'Ixxxx 

-155 

-2,092 

-263 

1,961 

3,205 

3,300 

95 

194 

2 

111 

2,008 

2,800 

^zzzz 

-6 

-13 

-23 

64 

11,303 

11,900 

81 

25 

-1 

43 

1,680 

1,600 

43 

82 

9 

17 

2,344 

2,500 

'Ixxzz 

98 

304 

40 

231 

3,294 

3,100 

Ji _ 

76 

-218 

-37 

-241 

6,230 

6,500 

^dcSHG 

^11 

85 

-344 

-42 

-337 

9,251 

9,900 

9,029 

±203 

96 

-811 

-49 

-538 

15,836 

17,500 
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Notice  that  a  minimum  ST0-3G  basis  underestimates  the  HF-SCF  dipole  polar¬ 
izability  by  over  a  factor  of  5,  with  double  zeta  (DZ)  still  being  in  error  by  ~30%. 
The  essential  role  of  polarization  functions  is  emphasized  by  the  6-31G+PD,  mean¬ 
ing  additional  diffuse  p  and  d  functions  on  C  and  H  {34),  and  the  (5s3p2d/3s2pld) 
POL+  basis  {14). 

The  basis  set  effect  is  amplified  dramatically  for  7,  where  four  products  of 
X,  y,  and  z  are  evaluated.  Note  that  in  inferior  basis  sets,  including  STO-3G  and  DZ, 
7|l ,  even  has  the  wrong  sign,  differing  from  the  better  converged  value  by  at  least  four 
orders  of  magnitude!  However,  once  a  few  polarization  functions  are  included,  even 
as  in  the  modest  6-31G+PD  basis,  eonvergenee  toward  the  Hartree-Fock  solution  is 
relatively  good.  The  remaining  differences  between  7|  and  7^^8110  which 

means  dc-induced  Second  Harmonic  Generation  and  Third  Harmonic  Generation, 
the  potential  experimentally  observed  quantities,  lie  in  the  frequency  dependence 
neglected  in  the  statie  HF  calculation.  We  will  discuss  that  aspect  later.  Obviously, 
it,  too,  is  numerically  important  in  providing  reliable  theoretical  predictions. 


Semi-empirical  Methods 

The  other  numerical  values  use  the  BMDO  (55)  and  spectroscopically  parameterized 
INDO/S  methods  (55),  (the  latter  is  also  known  as  ZINDO  (57)).  Like  nearly  all 
semi-empirical  methods,  INDO  assumes  an  underlying  minimum  Slater-iype-Orbital 
(STO)  basis  SCF  description,  which  is  close  to  that  of  STO-3G.  In  the  INDO  case, 
the  parameters  are  chosen  to  best  reproduce  the  minimum  basis  SCF  results.  Hence, 
if  this  were  done  sueeessfully,  results  about  on  the  level  of  STO-3G  would  be 
obtained,  clearly  a  level  far  inferior  to  that  required  for  hyperpolarizabilities,  and 
this  is  illustrated  by  the  observed  INDO  results. 

INDO/S  would  appear  to  have  a  little  better  chance  at  obtaining  reasonable  values. 
Despite  the  minimum  basis  set  description,  which  is  clearly  suspect,  the  spectroscopie 
parameterization  is  chosen  to  try  to  describe  excitations  to  the  low-lying  excited  states 
of  molecules  and  their  transition  moments  within  a  single  excitation  configuration 
interaction  (CIS)  description.  After  fitting  parameters  to  known  spectra  of  similar 
molecules,  the  method  is  expected  to  describe  related  molecules  reliably.  If  this  were 
true  for  all  excited  states,  from  the  SOS  formula  of  equation  31b,  obviously  the 
dipole  polarizability  would  be  well  described  as  it  requires  only  those  two  pieces  of 
information.  However,  the  low-lying  states  are  only  a  few  of  those  that  contribute 
to  the  polarizability.  Furthermore,  as  all  terms  have  the  same  sign,  even  though 
those  neglected  might  have  a  comparatively  high  exeitation  energy,  their  sum  total  is 
significant.  That  is  why  there  is  a  large  error  in  a,  particularly  the  component, 
compared  to  the  good  ab  initio  results.  Notiee  that  INDO/S  is  competitive  with  the 
ab  initio  DZ  description,  but  unlike  ab  initio  methods,  semi-empirical  results  eannot 
be  systematically  improved.  You  get  what  you  get! 

Once  you  expect  to  use  INDO/S  for  ^  or  7,  you  now  not  only  require  transition 
moments  between  ground  and  excited  states  but  between  excited  states  themselves. 
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Such  experimental  information  is  very  seldom  available,  much  less  for  the  plethora 
of  possible  excited  states,  to  offer  any  help  with  parameterization. 

Table  El  demonstrates  the  dramatic  failure  of  7y  for  ethylene,  where  the  sign 
is  wrong  and  it  is  off  by  4  to  5  orders  of  magnitude.  The  same  failure  happens 
for  For  the  small  molecules,  INDO/S  and  E^O  values  of  have  the  wrong 
sign  for  CO  and  NH3!  (See  also  (22)  for  other  examples.)  Consequently,  consid¬ 
ering  the  basis  set  limitations  inherent  to  semi-empirical  methods,  and  their  general 
inability  to  describe  transitions  between  excited  states,  only  generalizations  beyond  a 
minimum  basis  description,  and  further  and  more  severe  parameterization  explicitly 
for  hyperpolarizabilities,  should  enable  such  methods  to  offer  any  kind  of  quanti¬ 
tatively  reliable  results.  For  design  purposes,  the  hope  is  to  at  least  reproduce  the 
correct  trends  among  similar  molecules,  and  typically  for  only  one  dominant  axial 
component  rather  than  the  whole  tensor;  a  simpler  problem. 

The  other  widely  used  semi-empirical  methods,  like  those  in  MOPAC  (38,39) 
(namely  MNDO,  MINDO,  AMI,  PM3),  share  with  E4DO  a  minimum  basis  descrip¬ 
tion,  but  parameterization  is  attempted  to  be  made  directly  to  experimental  results 
like  dissociation  energies  and  other  properties,  instead  of  the  minimum  basis  SCF 
values.  In  particular,  unlike  INDO/S,  such  approaches  have  not  been  developed  to 
apply  for  excited  states  as  recommended  by  the  SOS  interpretation.  But  even  from  the 
energy  derivative  viewpoint,  which  pertains  to  INDO  without  the  “S”  as  well,  there  is 
still  little  reason  to  believe  that  such  methods  have  much  hope  of  reliably  describing 
qualities  as  sensitive  to  nuances  of  charge  distributions  as  are  hyperpolarizabilities. 
Several  MOPAC  examples  that  demonstrate  failures  are  presented  in  ref.  (22). 


Electron  Correlation 

up  to  now,  we  have  only  considered  HF  level  methods,  ab  initio  or  semi-empirical. 
After  ensuring  adequate  basis  sets,  there  are  two  particular  corrections  we  need  to 
consider:  one  is  the  frequency  dependence  (discussed  in  the  next  section)  and  the 
other  electron  correlation.  Both  would  be  essential,  including  their  mutual  coupling, 
to  offer  the  definitive  theoretical  study  for  the  purely  electronic  part  of  molecular 
hyperpolarizabilities.  To  initially  isolate  the  effects,  we  will  start  with  correlation 
corrections  for  static  hyperpolarizabilities. 

As  discussed  above,  HF  theory  makes  the  approximation  that  one  electron  moves 
in  an  average  field  of  n-1  other  electrons,  to  enable  replacing  the  two-particle  operator 
in  the  Hamiltonian  by  the  one-particle  Hamiltonian.  This  ignores  that 

i 

electrons  are  charged  species  causing  their  motions  to  be  instantaneously  correlated. 
Clearly,  the  correlation  of  electrons  bestows  an  additional  degree  of  stability  to  the 
molecule,  as  the  electrons  are  allowed  to  avoid  each  other,  and  this  effect  significantly 
contributes  to  the  molecule’s  charge  distribution  and  excited  states  description. 

In  equation  31,  we  derived  formulas  for  hyperpolarizabilities  based  upon  knowing 
the  exact  solutions  to  Ho,  which  is  the  Hamiltonian  in  the  absence  of  the  electric  field 
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perturbation,  £oW.  Such  exact  solutions  properly  include  all  two-electron  effects  in 
Tio,  meaning  their  eigenvalues  and  vectors  are  correlated.  This  should  be  contrasted 
with  replacing  Ho  by  the  Ho  =  Y^  f{i)  operator  as  employed  in  the  HF  theory.  These 

i 

additional  effects  of  correlation,  whether  used  from  the  energy  derivative  viewpoint 
or  the  SOS,  can  have  a  dramatic  effect  on  the  observed  results. 

We  can  consider  various  approaches  to  electron  correlation,  but  those  most  fre¬ 
quently  applied  to  hyperpolarizabilities  are  many-body  perturbation  theory  (MBPT) 
(4,5,40),  known  in  some  programs  as  MP)  and  coupled-cluster  (CC)  theory 
(7-9,14,41,42).  In  MBPT,  in  the  absence  of  the  electric  field,  electron  correla¬ 
tion  corrections  could  be  introduced  with  the  same  equations  derived  in  equation 
6-11,  where  Ho  is  chosen  to  be  Ho  =  ^d  the  correlation  perturbation 

y  =  replaces  W.  Then  =  ($„|y|V)(i)). 

This  defines  MBPT(2).  MBPT(2)  is  the  simplest  correlated  method,  consisting  of  the 
initial  contribution  due  to  double  excitations  from  occupied  to  unoccupied  orbitals 
(i.e.  ^fj).  MBPT  is  usually  applied  in  a  given  finite  order  n  [MBPT(n)],  where  n=4 
is  the  highest  frequently  used. 

In  the  absence  of  an  electric  field,  proper  treatment  of  correlation  in  hyperpolar¬ 
izabilities  requires  a  double  perturbation  approach  (43)  where  all  couplings  between 
V  and  £oW  are  allowed  (44),  with  £oW  applied  in  a  given  order  to  describe  the 
particular  polarizability,  and  V  is  preferably  included  in  all  orders.  A  straightforward 
double  perturbation  approach  is  possible,  but  usually  in  practice,  the  requisite  cou¬ 
pling  between  V  and  W  is  handled  in  two  other  ways.  The  first  way  is  by  correlating 
states,  and  then  adding  the  W  perturbation;  the  route  taken  in  equations  28  to  31:  This 
will  be  the  EOM-CC  route  described  in  Section  8.  Second,  we  can  take  the  viewpoint 
that  we  will  first  solve  the  HF  problem,  and  then  add  correlation  to  that  HF  solution 
straightforwardly,  except  that  Ho,  V,  $o  and  are  all  dependent  on  £o.  That 
is,  we  evaluate  correlation  corrections  to  the  electric  field  perturbed  Hamiltonian  as 
fe  =  fo  +  £oW  +  as  in  the  CHF  method,  so 

Ho{£o)  =  /e(^)  =  (jo  +  (*) 

i  i 

V  =  H(£o)  -  Ho{£o)  (36) 

E(£o)  =  EcEF  +  E^‘^\£o)  +  E^^\£o)  -h  •  •  • 

^CHf(£o)  =  ($o(£o)|^(£<,)|$o(£o)) 

and  all  evaluation  of  correlation  corrections  pertain  to  y(eo);  such  as  E^^\£o)  = 
($o(eo)|y(£o)|'0^^H^o)),  and,  similarly,  in  higher  orders  of  MBPT  (2,3,7,40).  Just  as 
in  HF  theory,  we  have  the  option  of  doing  this  analytically  or  as  a  finite  field. 

CC  theory  offers  a  natural  generalization  of  MBPT  that  sums  categories  of 
excitations  to  infinite  order.  For  example,  CCSD  means  all  single  and  double 
excitations  (i.e.  ’J'CCSD  =  exp  (Ti -f  T2)$o  where  $o  is  the  independent  particle 
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model  and  T\  and  T2  are  the  single  and  double  excitation  operators, 

=  (37) 

i<3 

a<Zb 

where  occupied  orbitals  i  and  j  are  replaced  by  unoccupied  orbitals  a  and  b, 

=  (38) 

i 

By  virtue  of  the  exponential  expansion,  'J'ccSD  also  contains  triple  excitations, 
T\T2^o,  quadruple  excitations,  T^^o,  etc.,  ensuring  a  highly  correlated  solution. 
CCSD  is  correct  through  MBPT(3)  with  many  additional  higher-order  effects  {9,28). 
The  CCSD(T)  model  includes,  in  addition,  a  non-iterative  {28)  evaluation  of  “con¬ 
nected”  triples,  T3$o  =  model,  correct  through  MBPT(4), 

i<i<k  •’  ■' 

a<6<c 

is  demonstrably  close  to  the  basis  set  limit  or  full  Cl  solution — ^the  ultimate  re¬ 
sult  {28).  We  can  simply  evaluate  the  CC  results,  just  as  described  above  for 
MBPT,  by  letting  V’CC  to  be  eo  dependent,  [(■0cc  =  exp  (T(eo))$o(f^o)).  along  with 
niso)  =  Ho  +  V{eo)  +  SoW  and  Bcci^o)  =  ($o(£o)|7^(£o)e^("»)|$4eo)]. 

Even  better,  because  of  the  advantage  of  analytical  evaluation  even  for  the  induced 
dipole,  we  might  describe  how  that  is  done  in  CC  theory.  MBPT  follows  as  a  special 
case  {J8).  Because  V’cc(f^o)  =  exp  (T(£o))^o(£o)  is  an  infinite  series,  we  do  not 
evaluate  the  induced  dipole  /?(£<,),  as 

m(£o)  =  (V’cc(£o)|W^|V’cc(eo))/(V’cc(£o)|V’cc(eo))  (39) 


because  the  expression  would  have  to  be  truncated.  Instead,  we  can  derive  the  form 
{18-28), 

ii{eo)  =  ($o(£o)|(l  +  A(eo))e-^("“)W^e-^("«’)|$o(£<,)) 


where  A  is  a  de-excitation  operator,  complementary  to  the  excitation  operator  T  and 
7pg  is  the  element  of  the  one-particle  “relaxed  density”  matrix  {28).  Both  A  and  T 
have  to  be  determined  from  the  CC  equations.  Subsequent  finite-field  differentiation 
relative  to  So  provides  the  various  polarizabilities  like  =  a,  =  2!/5,  etc. 

As  mentioned,  this  greatly  improves  precision  and  diminishes  the  number  of  distinct 
So  values  required. 


Table  IV  shows  the  effect  of  correlation  on  static  0, 0)  and  Xy^^O;  0, 0, 0). 

(2') 

Note  the  factor  of  2  change  for  Xy  of  NH3,  and  the  order  of  magnitude  and  sign 

change  for  H2S.  Similarly,  x||^^  changes  substantially.  This  large  correlation  effect  has 
been  observed  since  the  initial  correlated  studies  of  molecular  hyperpolarizabilities 
(2,5),  and  makes  it  apparent  that  “predictive”  ab  initio  methods  for  hyperpolarizabil¬ 
ities  must  include  electron  correlation. 
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Table  IV  Comparison  of  SCF  and  Correlated  Static  Hyperpolarizabilities.^^ 


“  =  (Aii +/^j>)/5  with  summation  convention.  =  /5||/2,  so 

conversion  to  a.u.  is  4.3195  x  10"^^  esu/molecule/a.u. 

Ty  =  {liijj  +7ijij  +7wi)/15  with  summation  convention.  Xy^^  =  7y/6,  so 
conversion  to  a.u.  is  8.395  x  esu/molecule/a.u. 

®  Basic  sets  listed  in  Tables  Vn  and  Vin. 

Frequency  Dependence 

All  experiments  are  frequency  dependent,  and  frequency  dependence  introduces  many 
different  processes  that  become  the  same  in  the  static  limit. 

Instead  of  the  static  field  perturbation,  Eo  ■  W,  consider  the  expansion  of  the 
induced  dipole  analogous  to  that  in  equation  24,  for  a  time-dependent,  oscillatory 
field,  £  =  So  +  cos  ut, 
fii  =  (v?(e,f)|rj|\t'(£,i)) 

=  f^iio)  -f  aij(o,o)£oj  +  uj)eoj  cos  wt 

d"  ~^Pijk{p\  ^)^uij^u>k  (41) 

+  COS  Ujt  +  -Pijfi[  —  2uj-,(V^u)£uj£i_jk  COS  ut 

+  ... 

Now,  in  addition  to  the  static  terms  we  have  previously  considered,  we  obtain  a 
number  of  terms  that  correspond  to  different  incoming  and  outgoing  frequencies.  For 
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example,  2a;;a;,ti;)  corresponds  to  second  harmonic  generation  (SHG)  with 

incoming  frequencies  a;  +  w  resulting  in  an  outgoing  frequency,  2a;.  Similarly, 
corresponds  to  the  dynamic  polarizability,  and  Pijk{o-,<jO,-Lj)  is  called 
optical  rectification  (OR),  and  u,  o)  corresponds  to  the  electro  optic  Pockels 

effect  (EOPE).  If  we  allowed  the  frequencies  to  be  different,  /3(— a;i  —  uJ2‘,uji,i02) 
would  correspond  to  two-wave  mixing.  Similarly,  many  such  components  occur  in 
7.  A  summary  of  them  is  shown  in  Table  V.  All  become  the  same  in  the  static  limit; 
hence,  without  inclusion  of  frequency  dependence  in  the  quantum  mechanical  method, 
we  cannot  distinguish  between  the  different  processes.  Obviously,  we  can  obtain  each 
of  these  quantities  from  an  appropriate  derivative. 


d^Pije) 


=  ^Aifc(-o;a;,-a;) 


(42) 


etc.  Just  as  in  the  static  case,  that  derivative  can  be  further  related  to  a  quasi-energy 
derivatives  (45,46)  from  the  Frenkel  variational  principle  (11). 


Table  V  Representative  Non-Linear  Optical  Processes,  with  Corresponding  Resultant 
(uja)  and  Incident  (a;i,a;2  ...  )  Frequencies. 


Non-Linear  Optical  Process 

-wi 

-ajj 

First  Static  Hyperpolarizability 

0 

0 

- 

Second  Harmonic  Generation  (SHG) 

-2ct^ 

UJ 

(Jj 

- 

Electro-Optics  Pockels  Effect  (EOPE) 

~0J 

0 

(Jj 

- 

Optical  Rectification  (OR) 

0 

UJ 

-UJ 

- 

Two- Wave  Mixing 

-(^^l+^z) 

Wx 

Second  Static  Hyperpolarizability 

0 

0 

0 

0 

Third  Harmonic  Generation  (THG) 

-3a; 

(jJ 

(jLf 

Intensity-Dependent  Refactive  Index 
(IDRI) 

-UJ 

to 

(jJ 

Optical  Kerr  Effect  (OKE) 

-UJl 

Ui 

(4J2 

’OJ2 

D.C.-Induced  Optical  Rectification 
(DCOR) 

0 

0 

U 

-U 

D.C.-Induced  SHG  (DC-SHG) 

-2a; 

0 

(jJ 

UJ 

Electro-Optic  Kerr  Effect  (EOKE 

-ijJ 

U) 

0 

0 

Three- Wave  Mixing 

-U)(J 

a;i 

U)2 

D.C.-Induced  Two-Wave  Mixing 

-(a;i-i-a;2) 

0 

a;i 

U)2 

The  first  way  to  augment  a  static  calculation  with  some  measure  of  frequency 
dependence  is  to  recognize  that  it  may  be  rigorously  shown  (47)  that  for  low 


frequencies, 


=  /3(o;  o,  o)  (l  4-  AJi  +  Bu}\  ^ - )  (43) 

where  ix)l  =  (^^+t^i+co2  A  and  B  are  unknown  constants.  Without  an  independent 
evaluation  of  A  and  B,  the  most  we  can  conclude  are  the  ratios  of  the  various  dispersion 
effects. 

For  example,  neglecting  the  smaller  quartic  term,  the  SHG  and  HOPE  values  are 

/3(-2uj;aj,u>)  =^(o;o,o)(l  +  6Au^) 
w,  o)  —  p[o]  o,  o)  (l  +  2Au>^'j 


showing  that  their  ratio. 


P(—2u;;  ojyUj) 


1  +  6Auj^ 

1  +  2Aw2 


~  1  +  4Aa;^  —  SA^ui'^  H - 


(45) 


Obviously,  SHG  has  a  much  greater  dispersion  effect  than  EOPE.  Similarly,  we  have 
for  7, 


=  7(0;0,0,0)(l  +  A'a;|  +  S'w^  +  ---)  (46) 

The  components  of  7,  OKE,  IDRI,  DCSHG  and  THG  correspond  to  the  values 
1  plus  2A'uj^,  4A'u}^,  GA'uj^  and  12A'a;^  ,  respectively,  providing  similar  ratios. 
Obviously,  the  degree  of  dispersion  follows  the  order  THG  >  DCSHG  >  IDRI  > 
OKE.  However,  we  must  at  least  know  the  constants  to  relate  the  various  processes, 
particularly,  to  static  quantities.  It  appears  the  only  way  at  present  to  obtain  any 
quantitative  relationship  is  to  evaluate  the  frequency  dependent  quantum  chemical 
results. 

The  first  such  viable  approach  is  the  time-dependent  Hartree-Fock  (TDHF) 
method.  Just  as  HF  theory  provides  the  simplest  approach  for  static  polarizabilities, 
TDHF  provides  the  simplest  ab  initio  solution  to  the  time-dependent  Schrddinger 
equation  =  and,  thereby,  frequency  dependent  processes.  Just  as  in 

the  static  case,  we  can  impose  a  single  determinant  approximation  for  ^{t)  «  $o(t), 
and  by  applying  the  time-dependent  (Frenkel)  variational/principle  obtain  the  TDHF 
(or  RPA)  equations.  The  REA.  equations  date  from  long  ago  (48);  however,  this 
only  pertains  to  excitation  energies  and  a.  In  1986,  we  made  the  generalization  to 

7)  i)  and  showed  that  analytically  we  can  evaluate  the  frequency  de¬ 
pendent  hyperpolarizabilities  from  the  derivatives  in  equation  36  in  any  order  (11) 
giving  ready  TDHF  access  to  hyperpolarizabilities.  Others  have  since  implemented 
equivalent  TDHF  approaches  for  p  (45)  and  for  p  and  7  (32).  Recently,  the  unre¬ 
stricted  Hartree-Fock  (open  shell)  generalization  has  been  made  (49).  Figs.  1  and  2 
graphically  illustrate  the  dispersion  behavior  for  various  processes,  which  reflect  the 
relative  numerical  proportions. 
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ate 


Figure  1 
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Frequency  (a.u.) 


Frequency(a 


Table  VI  shows  the  effect  of  dispersion  in  HF-level  calculations  of  frequency 
dependent  polarizabilities  to  SHG  and  dcSHG  processes.  For  the  average  change 

is  20%  and  Xy  is  26%.  Obviously,  from  the  static  values  the  other  dynamic  processes 
can  show  greater  and  lesser  effects.  Just  as  usual,  we  can  take  the  derivative  or  SOS 
viewpoint  in  TDHF.  In  the  latter  case,  the  excited  TDHF  states  are  the  usual  RPA 
ones  (24). 

Table  VI  Comparison  of  Static  and  Dynamic®  Hartree  Fock  Hyperpolarizabilities.^ 


x[|^^  esu/molecule) 

Xn^^  esu/molecule) 

SCF 

TDHF 

SCF 

TDHF 

H2 

- 

46 

53.6 

N2 

- 

- 

61 

69.0 

CO2 

- 

67 

76.6 

C2H4 

- 

- 

546 

832 

CO 

9.1 

10.5 

85 

102 

HF 

-2.3 

-2.5 

27 

30 

HCl 

-1.3 

-1.6 

213 

270 

H2O 

-4.7 

-5.4 

85 

100 

NH3 

-6.5 

-9.3 

200 

280 

H2S 

0.58 

0.64 

470 

690 

®  Values  by  TDHF  for  SHG  and  dcSHG  at  694.3  nm. 
Basis  sets  listed  in  Tables  VII  and  Vin. 


At  this  point,  we  have  a  procedure  based  upon  TDHF  to  evaluate  the  dispersion, 
and  a  procedure  to  add  electron  correlation  to  static  hyperpolarizabilities.  Clearly, 
both  are  critical  in  obtaining  predictive  values.  Hence,  guided  by  the  fact  that  equa¬ 
tions  37  and  40  are  exact  in  the  low-frequency  limit,  it  makes  sense  to  use  a  percentage 
TDHF  dispersion  correction  (13)  to  augment  static,  correlated  hyperpolarizabilities; 
namely, 

Phf(0;  0, 0) 

=  7(0; 0,0,0)  X  TrDHF(^,; 

7hf(0;  0,0,0) 

By  equating  the  percentage  correction  to  Au)\  at  a  particular  frequency,  a  value  of  A 
could  be  extracted,  as  well;  or  fitting  to  several  different  processes,  A  and  B.  In  the  few 
cases  where  frequency  dependent  correlated  results  have  been  obtained  (20,50,51), 
the  percentage  TDHF  dispersion  estimate  has  been  well  supported.  However,  it  is 
clear  that  if  the  TDHF=RPA  result  for  the  excitation  energies  are  poor,  then  the  slope 
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of  the  curve  in  Figs.  1  and  2  will  have  to  eventually  change  to  be  able  to  approach 
the  different  asymptotic  values  of  the  excitation  energy. 

Using  this  decoupled  TDHF  dispersion,  static  correlation  procedure,  we  obtain  the 

results  shown  in  Tables  Vn  and  Vin.  All  fall  within  10%  error  of  the  experimental 

(2) 

result  except  for  FH,  whose  errors  are  28%  for  Xy  and  24%  for  Xy  ,  for  C2H4, 

(3) 

whose  Xy  is  too  large  by  13.4%  and  N2,  and  CO,  whose  error  is  a  modest  11%. 
We  have  discussed  FH  at  length  elsewhere  {7,14),  so  we  will  not  repeat  that  here, 
except  to  say  that  we  dispute  the  experimental  values,  expecting  an  error  in  its 
determination.  Extensive  studies  of  the  convergence  of  FH’s  hyperpolarizabilities 
to  basis  set,  correlation,  frequency  dependence  and  vibrational  corrections  {14),  show 
no  convergence  to  the  experimental  values.  Instead,  we  propose  that  the  correct  value 
for  Xy^^  is  -3.6  ±  0.3x10'^^  esu/molecule  and  Xy^^  is  55  +  5x10’^^  esu/molecule.  For 
the  similar  HCl  molecule,  our  CCSD(T)  results  fall  within  our  goal  of  a  10%  error. 

Table  VII  Theoretical  vs.  Experimental'’  Hyperpolarizabilities  of  Molecules. 


.w(2)  _  1  /3SHG 

^11  ~  2ni 

TDHF 

MBPT(2) 

CCSD 

CCSD(T) 

EXPERIMENT 

CO 

10.5 

11.2 

11.4 

11.7 

12.9  ±1.4 

HF 

-2.5 

-3.3 

-3.2 

-3.4 

-4.70  ±0.41 

HCl 

-1.6 

-4.0 

-3.3 

-3.8 

-4.22  ±0.50 

H2O 

-5.4 

-8.8 

-8.2 

-9.1 

-9.4  ±0.4 

NH3 

-9.3 

-20.1 

-18.4 

-21.2 

-20.9  ±0.5 

H2S 

0.6 

-5.0 

-3.4 

-4.5 

-4.3  ±0.9 

“Value  corrected  for  the  dispersion  effect  at  694.3  nm  using  the  formula. 


TDHF(a;) 
TDHF(£^  =  0) 


The  calculations  are  performed  with  basis  sets  [5s3p2d]  for  C,  N,  O  and  F;  [7s5p2d] 
for  S;  and  [3s2pld]  for  H.  Lone-pair  functions  are  added  for  HF,  H2O  and  NH3.  For 
HCl,  basis  is  [8s6p3dlf]  for  Cl  and  [3s2pld]  for  H.  All  molecules  at  experimental 
geometries  and  there  is  no  estimate  of  vibrational  corrections.  H2S  values  in  lone-pair 
augmented  basis  are  1.0,  -4.4,  -2.8  and  -3.8,  respectively  {14).  The  lone-pair  basis 
HCl  values  are  -1.0,  -3.1,  -2.6  and  -3.0. 

'’Value  obtained  by  dc-induced  Second  Harmonic  Generation. 
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Table  VIII  Theoretical®  vs.  Experimental'^  Hyperpolarizabilities  of  Molecules. 


v(3)  _  UdcSHG 

_ _ _ _ _ X|[ . z.e^ll _ 

TDHF 

MBPT(2) 

CCSD 

CCSD(T) 

EXP 

H2 

54 

58.2 

59.3 

59.3 

60.5^ 

N2 

69 

88.7 

90.6 

96 

86.6  ±1.0 

CO2 

77 

111.5 

107.9 

110 

111.9+1.3 

C2H4 

822 

960 

820 

860 

758  ±17 

CO 

102 

151 

149 

160 

144  ±4 

HF 

30 

52 

49 

53 

70  ±10 

HCl 

270 

364 

352 

374 

347  ±15 

H2O 

100 

180 

170 

180 

194  ±10 

NH3 

280 

460 

430 

470 

511  ±9 

H2S 

690 

910 

870 

930 

865  ±22 

® Value  corrected  for  the  dispersion  effect  at  694.3  nm  using  the  formula. 


TDHF(a;) 
TDHF(a;  =  0) 


The  calculations  are  performed  with  basis  sets  [5s3p2d]  for  C,  N,  O  and  F;  [7s5p2d] 
for  S;  and  [3s2pld]  for  H.  The  lone-pair  functions  are  added  for  HF,  H2O  and 
NH3.  For  HCl,  basis  is  [8s6p3dlf]  for  Cl  and  [3s2pld]  for  H.  All  molecules  at 
experimental  geometries  and  there  is  no  estimate  of  vibrational  corrections.  The 
lone-pair  augmented  H2S  basis  values  are  731,  970,  930  and  980,  respectively.  The 
lone-pair  basis  HCl  values  are  274,  369,  356  and  378. 

'’Value  obtained  by  dc-induced  Second  Harmonic  Generation. 

‘’Exact  electronic  value  [D.  M.  Bishop,  J.  Pipin  and  S.  M.  Cybulski,  Phys.  Rev. 
A43,  4845  (1991)].  The  experimental  value  of  67.2  includes  a  significant  vibrational 
effect. 


For  C2H4  and  N2,  we  can  consider  more  recent  or  rescaled  experimental  values 
(22).  These  are  76.6  ±  17  and  88.8  ±  1,  respectively.  This  does  reduce  the  C2H4  error 
to  12.3%  and  that  for  N2  to  7.9%.  Using  the  experimental  results  in  the  tables,  the 
average  error  for  CCSD(T)  is  9.4%.  Excluding  FH,  it  becomes  7.3%.  Similarly, 

the  average  error  for  is  9.4%,  or  5.7%  without  the  FH  example,  falling  within 
our  10%  error  target  objective.  The  various  contributions  are  shown  in  Table  DC. 
Clearly,  correlation  and  frequency  dependence  are  critical,  and  CCSD(T)  is  better 
than  CCSD.  However,  it  is  gratifying  that  MBPT(2),  which  is  a  lot  cheaper  than 
CCSD(T),  maintains  about  a  10%  error.  This  level  of  correlation  can  be  applied  to 
much  larger  molecules  (52)  than  can  CCSD(T). 
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Table  IX  Percent  Error  of  Hyperpolarizabilities  at  Various  Levels  Compared  to 
Experiment. 


(2) 

X|| 

X|| 

HE  iuj=0) 

57 

46 

TDHF  (w) 

50 

31 

MBPT(2) 

12  (9)* 

9  (7)* 

CCSD 

19  (16)* 

8  (6)* 

CCSD(T) 

9  (6)* 

9  (7)* 

*Without  FH  example. 

Correlated  Frequency  Dependent  Polarizabilities 

Despite  the  success  of  the  decoupled  correlated/TDHF  results  shown  above,  the  most 
rigorous  method  would  include  the  full  coupling  between  correlation  and  dispersion. 
There  have  only  been  three  attempts  of  this  type  for  hyperpolarizabilities:  a  second- 
order  [MBPT(2)=MP(2)]  level  method  (50);  a  multi-configurational  (MCSCF)  linear 
response  approach  (57);  and  our  EOM-CC  method  (20).  The  first  two  are  devel¬ 
oped  from  the  energy  (or  quasi-energy)  derivative  viewpoint,  while  the  latter  refers, 
conceptually,  to  the  SOS  expressions  in  equation  31.  In  other  words,  EOM-CC  pro¬ 
vides  excited  states,  their  excitation  energies,  and  generalized 

transition  moments,  from  which  the  SOS  expressions  could  be  formally 

constmcted.  That  form  is  particularly  convenient,  since  frequency  dependence  can 
be  trivially  added  to  such  an  SOS  expression.  For  example,  for  some  frequency  u). 


Wfc  +  w 


U>k  —  (jlJ 


] 


(50) 


When  oj  =■  Eo  -  we  have  a  pole,  whose  residue  is  the  dipole  strength 

The  basic  idea  of  EOM  is  very  simple  (53).  Consider  the  solutions  to  the 
Schrbdinger  equation  for  an  excited  state,  0^  and  for  the  ground  state,  'ipo, 

Hoi>o  =  Eo-ipo  (51a) 

Wo#  =  (5Ib) 


Now  we  will  choose  to  write  ,  where  Ilk  is  an  operator  that  creates 

excitations  from  If  we  limit  ourselves  to  single  and  double  excitations,  we 


28 


will  define  an  EOM-CCSD  approximation.  Inserting  the  excited  state  expression  into 
equation  49b,  multiplying  equation  49a  by  from  the  left  and  subtracting,  we  obtain 


ilioTZk  -  T^k'^o)i’o  =  -  Eo^TZki’o 

[7^0;  =  uJk'E.k'fpo 


(52a) 

(52b) 


with  the  commutator  between  Ho  and  Tlk,  [i.e.  expression  gives  the 

“equation-of-motion,”  from  the  obvious  connection  with  Heisenberg’s  form.  To 
introduce  CC  theory,  we  simply  choose  for  the  correlated  ground  state,  V'o  =  "^cc  = 
exp(ri  +T2)#o.  As  TZk  and  !„  are  all  excitation  operators,  [7^fc,T„]  =  0,  and  we 
can  commute  the  operators  to  give, 

[e~^Hoe'^,TZk]^o  =  [Ho,'R,k]^o  —  (^kE^k^o  (53) 


recognizing  that 


=  <•.  +  X;  n-*?  +  E  '■“/ 

i.a 

a<b 

the  coefficients  {rf,rfj}  are  to  be  determined  by  the  matrix  equation, 

Htfe  =  TkLOk 


(54) 


(55) 


Since  Ho  is  not  Hermitian,  we  also  have  left-side  eigenvectors,  which  unlike 

Ek^o,  correspond  to  de-excitation  processes.  The  Ck  and  TZi  states  are  biorthog- 
onal,  {^o\EkEi\^o)  =  ^ki-  Their  matrix  equation  is,  l^Ho  =  tokh-  The  corre¬ 
sponding  generalized  transition  moments  are  obtained  from  = 

{^o\EkfiEi\^o)-  Notice  ($o|An7^/|^o)  7^  {^o\EifiEk\^o)  since  the  operator  fi  is 
not  Hermitian.  The  observable  quantity,  however,  is  the  dipole  strength  (the  prod¬ 
uct),  not  the  transition  moment  itself. 

Now  if  we  return  to  equation  48,  in  terms  of  EOM-CC  solutions,  we  have 


(56) 


-f- 


Wfe  -  a> 


Notice,  we  used  the  fact  that  (#o|£o  =  (^o|(l  +  A),  which  is  the  same  A  operator 
introduced  in  equation  40,  and  that  =  1- 

Though  informative  in  this  form,  we  would  still  have  to  truncate  the  SOS. 
To  avoid  any  such  truncation,  we  need  to  recognize  that  all  EOM-CC  states  are 
ultimately  represented  in  terms  of  their  underlying  single  and  double  excitations,  i.e. 
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the  expression  in  equation  54.  If  we  collectively  represent  all  as  |h),  it 

may  be  shown  (79)  that  equation  48  may  be  written  as 

=  ((#„|(1  +  A)(fi  - 
+(4,1(1 +  A)(ri-(0»ri‘>|4„)) 

where  (fi)  is  the  generalized  expectation  value  like  that  in  equation  40,  without 
£o  dependence).  Analogous  to  ordinary  perturbation  theory  [equations  28  to  30], 
the  first-order  perturbed  correlated  wavefunction  is  given  in  terms  of  the  resolvent 
operator  matrix  Ro, 

{h|T^'*|4„)  =  (h|Bcc  ±1^  -  H„|h)-"(h|fi|4„)  =  R„(h|f,  |4,)  (58) 


In  practice,  instead  of  inversion,  we  solve  the  very  large  linear  equation, 

(h|Ecc  io;  -  ?7,|h)(h|Ti^’^|$o)  =  (h|fj|$,)  (59) 

at  a  given  value  of  oj  for  the  coefficients.  Hence,  we  can  evaluate  the  SOS 
dynamic  polarizability  without  any  truncation  {19,20)\ 

From  the  above  evaluation  of  u),  we  can  obtain  the  hyperpolarizabilities 

for  the  optical  Kerr  effect  (OKE)  as  follows.  ^{—oj\uj,o),  sometimes  also  called  the 
HOPE,  is  the  second-order  hyperpolarizability  obtained  from  a  Kerr  effect  experiment, 
while  the  EOKE  corresponds  to  7(-a;;  w,o,o).  Because  of  the  static  fields  in  both 
processes,  we  are  able  to  obtain  P{-<jj-,u;,o)  and  7(-a;;a;,o,  o)  from  finite-field 
differentiation  of  a{—u};(jj).  That  is,  we  evaluate  aij{—uj-,uj)  analytically,  using 
equations  56  and  57,  where  we  use  the  perturbation  £„  •  f  instead  of  just  f.  Then  we 
obtain  a(—(jt};ujo,£o),  from  which 


l3^{-u\uj,o)  = 
7*(-a;;w,o,o) 


da{—u}\uj\£o) 

dso 

del 


(60a) 

(60b) 


Results  are  shown  in  Table  X  for  NH3  and  trans-butadiene  as  a  function  of  frequency, 
for  EOM-CC  and  TDHF.  The  usual  very  large  effect  of  correlation  accounts  for  the 
much  larger  magnitude  for  the  EOM-CC  values,  while  the  comparative  dispersion 
values  are  indicated  as  the  percent  dispersion  in  parentheses.  For  NH3,  there  is 
about  10%  greater  dispersion  as  measured  by  EOM-CC  compared  to  TDHF  at  the 
high  frequencies  (0.1  a.u.),  but  not  at  the  low  frequency  (0.0656  a.u.)  value  used 
in  the  dcSHG  experiments  we  previously  described.  For  butadiene,  the  percentage 
dispersion  is  close,  but  the  slightly  smaller  value  helps  to  reduce  our  calculated  values 
to  be  in  somewhat  better  agreement  with  experiment.  Unlike  most  other  molecules 
we  have  studied,  TDHF  results  for  ethylene  and  butadiene  are  fortuitously  close  to 
experiment,  while  correlation  hurts  the  agreement.  The  origin  of  this  is  not  yet  clear. 
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but  the  fact  that  the  correlated  dispersion  is  smaller  than  that  for  TDHF  could  be 
ascribed  to  the  TDHF=RPA  excitation  energies  being  too  low,  causing  the  curve  in 
Fig.  2  to  rise  too  quickly  to  approach  the  wrong  asymptotic  values.  For  multi-bonded 
molecules  like  ethylene  and  butadiene,  the  restricted  HF  result  is  not  triplet  stable, 
although  this  only  prohibits  RPA  from  correctly  describing  triplet  excited  states,  the 
RPA  singlet  excitations  tend  to  be  lower  than  experiment;  contrary  to  that  for  most 
systems.  This  may  partially  account  for  TDHF  giving  higher  percentage  dispersion 
corrections. 

Table  X  Kerr  Effect  Tensors  for  NH3  and  C4H6  Calculated  at  Different  Frequencies 
(in  au).“‘ 


Frequency 

(au/nm) 

0 

0.043/1060 

0.0656/694.3 

0.1/455.6 

NHs^  (POL+) 

pkb 

-35.9 

-14.7 

“ 

-41.5  (15.7%) 
-16.6  (12.7%) 

-51.5  (43.5%) 
-19.7  (33.8%) 

ryk  C 

4136.63 

2404.99 

- 

4703.00  (13.7%) 
2646.36  (10.0%) 

5711.56 

(38.1%) 

3039.95 

(26.4%) 

C4H/  (631G+PD) 

lxxxx{-<^',<^,0fi) 

44000  (6.9%) 
25571  (8.7%) 

48100  (16.8%) 
28733  (22.2%) 

- 

jk 

Ml 

21900  (5.7%) 
15794  (6.6%) 

23700  (14.1%) 
17277  (16.7%) 

- 

®  The  numbers  in  upper  and  lower  rows  are  evaluated  by  EOM-CC  and  TDHF, 
respectively,  x  is  the  Cj  molecular  axis, 
b  /3^  = 

^  ~  liijj) 

The  numbers  in  upper  and  lower  rows  are  evaluated  by  EOM-CC  and  TDHF, 
respectively.  The  x  component  corresponds  to  the  longitudinal  molecular  axis  and 
the  z  component  is  perpendicular  to  the  molecular  plane. 

We  have  not  discussed  vibrational  contributions  to  predictive  studies  of  hyper¬ 
polarizabilities,  but  these  can  sometimes  be  important  {21).  In  equation  48,  e.g.,  we 
could  have  contributions  from  all  vibronic  states  indexed  by  k,  instead  of  just  the 
electronic  ones.  For  an  optical  frequency,  uj,  which  is  much  greater  than  a  vibrational 
energy,  such  slight  changes  in  would  have  negligible  numerical  value.  How¬ 
ever,  for  NLO  processes  that  involve  static  fields,  the  SOS  formulas  will  have  some 
denominators  without  an  additional  large  u  value,  causing  the  vibronic  changes  in 
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to  be  more  significant  to  the  final  result.  The  vibrational  energy  levels  can  be 
substantially  perturbed  by  such  an  electric  field  as  we  have  shown  numerically  for  FH 
and  H2  (10).  Because  of  the  static  fields  in  OKE,  this  is  an  example  where  attention 
needs  to  be  paid  to  such  effects.  Assisted  by  a  determination  of  /S^and  (54),  we 

can  extract  from  our  calculations  predictions  for  and  7*  =  that 

could  be  compared  with  experiment.  For  NH3,  MBPT(2)  values  for  and  7  k  are 
3.8  and  135  a.u.  at  0.07  a.u.  This  suggests  that  «  —36,  7*  «  4800  at  ^=0.0656. 

For  butadiene,  there  is  a  larger  (SCF),  7^=1395  to  1762  a.u.,  [55]  giving  «25000. 

Future  Extensions 

There  are  a  couple  of  fairly  obvious  extensions  that  should  be  made  in  future 
theoretical  work  for  NLO  material  design.  In  the  short  term,  we  obviously  need 
to  generalize  analytical  frequency  dependent  EOM-CC  for  all  tbe  components  of 
^  and  7.  Also,  considering  the  good  accuracy  of  MBPT(2)  level  correlation,  purely 
analytical,  frequency  dependent  versions  for  P  and  7  are  strongly  recommended, 
that  should  remove  the  current  constraints  (50).  Other  routes  to  partitioned  EOM- 
CC  approximations  that  are  operationally  second-order,  can  be  envisioned  (56)  and 
should  be  pursued. 

Obviously,  it  would  also  be  nice  to  be  able  to  treat  hyperpolarizabilities  for 
molecules  in  solution.  Several  sucb  solvation  methods,  ranging  from  continuous 
reaction  fields  to  more  detailed  solvation  models  are  becoming  available  (57).  These 
should  help  in  sorting  out  the  large  discrepancies  among  results  from  solvation 
experiments  (EFISH)  (58). 

In  the  longer  term,  we  need  theoretical  methods  comparable  to  that  presented  for 
small  molecules  that  are  applicable  to  extended,  polymeric  systems  (59).  The  first 
such  approaches  should  employ  periodicity,  with  future  extensions  directed  at  the 
inclusion  of  impurity  effects. 

Today,  it  is  not  possible  to  use  analytical  gradient  techniques  with  correlation  to 
move  atoms  around  in  polymers,  as  it  is  for  molecules.  Nor,  are  there  the  quality 
ab  initio  methods  for  band  gaps,  and  excited  states,  and  polarizabilities  as  there  are 
for  molecules.  Clearly,  developing  tbe  tools  for  rational  NLO  polymer  design  should 
have  a  high  priority. 

As  current  high-level  ab  initio  methods  will  eventually  encounter  limitations, 
even  for  periodic  systems,  simplified  techniques  should  be  pursued,  simultaneously. 
The  questionable  reliability  of  semi-empirical  MO  theory  suggests  that  a  better  “semi- 
empirical”  approach  is  likely  to  be  offered  by  modern  density  functional  theory  (DFT). 
Although  DFT  has  a  rigorous  base,  in  application  it  is  semi-empirical.  Such  methods 
are  well  known  for  extended  systems,  and  have  decided  computational  advantages, 
compared  to  ab  initio  correlated  methods,  but  they  have  not  yet  been  demonstrated 
to  provide  comparable  results  to  those  presented  in  this  chapter.  In  fact,  one  paper 
says  that  Kohn-Sham  DFT  does  not  work  for  molecular  hyperpolarizabilities  (60).  We 
have  considered  other  DFT  variants,  however,  and  find  that  competitive  results  can  be 
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obtained  (61).  The  significant  computational  advantages  of  DFT  make  this  a  profitable 
area  for  study.  Frequency  dependent  approaches  need  to  be  developed,  however. 
Also,  the  conventional  wisdom  is  that  DFT  does  not  admit  treatments  for  excited 
states.  Exploiting  the  equivalent  derivative  viewpoint  should  avoid  any  such  formal 
restrictions  for  polarizabilities.  Also,  the  ultimate  limitation  of  applied  DFT  methods, 
like  semi-empirical  MO  methods,  is  that  there  is  no  way  to  systematically  converge 
to  the  exact  result.  New  methods  that  combine  elements  of  ab  initio  correlated  theory 
with  DFT  methods  will  be  forthcoming  and  might  alleviate  this  failing. 
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